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During the 2008 run RHIC provided high luminosity in both p+p and d+Au collisions at 
y/SNN = 200GeV. Electromagnetic calorimeter acceptance in STAR was enhanced by the 
new Forward Meson Spectrometer (FMS), and is now almost contiguous from — 1 < n < 4 
over the full azimuth. This large acceptance provides sensitivity to the gluon density in the 
nucleus down to x ~ 10~ 3 . Measurements of the azimuthal correlation between a forward 
7T° and an associated particle at large rapidity are sensitive to the low- a; gluon density. Data 
exhibit the qualitative features expected from gluon saturation. A comparison to calculations 
using the Color Glass Condensate (CGC) model is presented. 



1 Introduction 

Forward acceptance at STAR has been drastically improved in the last few years with the 
development and commissioning of the Forward Meson Spectrometer (FMS). The FMS is a high 
granularity electromagnetic calorimeter whose purpose is to measure photons from decays of 
forward neutral mesons. It extends the fully azimuthal, electromagnetic capability of STAR 
into the pseudorapidity region 2.5 < r\ < 4.0, making the overall coverage almost hermetic in 
the wide —1.0 < rj < 4.0 interval. This allows measurements of correlations of different species 
of forward- and mid-rapidity particles, over a broad Ary x Aip range. 

One of the main purposes of the FMS is to characterize the px dependence of the di-pion 
azimuthal correlations as a means to search for the onset of saturation effects^ in the nuclear 
gluon distribution. It is believed that, as the energy grows, non- linear effects must be included 
in the nuclear wave-function in order to tame the otherwise divergent rise of the gluon density. 
At very low values of the longitudinal momentum fraction x of the probed gluon, as accessed 
in high-energy collisions, the occupation numbers become large, allowing gluon recombination 
that eventually leads to saturation. Saturation is expected to be revealed at RHIC in d+Au 
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Figure 1: On the left: schematic view of the STAR hall. Gold beam is coming from the West (right in figure) 

side. On the right: uncorrected coincidence probability versus azimuthal angle difference between two forward 

neutral pions. Different pr cuts are indicated inside the figure. The left (right) column shows p+p (d+Au) data 

with statistical errors. Data are fit with a constant plus two gaussian functions (in red). 

collisions, where large densities of nuclear gluons are probed with a much simpler final state than 
heavy ion collisions. The FMS, facing the deuteron beam direction, provides measurements of 
forward hadron production in the rapidity region that selects small- x gluons in the nucleus. It 
probes x down to x ss 10~ 3 for inclusive particle production at (rj) ss 3 at -v/saF/v = 200 GeV^, 
well into the range where saturation effects are expected to set in. 



2 Azimuthal correlation analysis 

Models try to describe forward hadron production from dense nuclear targets by including non- 
linear effects. One approach^ extends perturbative QCD factorization by adding contributions 
of quarks scattering coherently off multiple nucleons, leading to an effective shift in the gluon 
x. In the Color Glass Condensate (CGC) model^^ the low-x gluon density is saturated and 
non-linearities are treated classically. In both cases the 2—7-2 picture of elastic parton scattering 
process is replaced by a 2 — > many picture, where the probe scatters coherently off the dense 
gluon field of the target, which recoils collectively. In the CGC model, the transverse momentum 
of a jet produced by the large momentum parton (most likely a valence quark) in the deuteron 
is balanced by many gluons in the nucleus. To quantify this effect, it is instructive to compare 
the azimuthal angular correlation (Ay?) between the forward tt° and coincident hadrons. In p+p 
collisions, one expects to see a peak centered at Atp = it, representing the elastic back-to-back 
scattering contribution. In forward d+Au collisions, on the other hand, non-linear contributions 
are expected to cause a loss of correlation between the two particles, leading to a broadening of 
the back-to-back peak, and eventually to its disappearance (monojet). 

A systematic plan of azimuthal correlation measurements has been pursued at STAR using 
forward prototype calorimeters (Forward Pion Detector, FPD/FPD++) and, more recently, the 
FMS. The objective has been to probe if the boundaries of the saturation region are accessible at 
RHIC energies^ and establish the effect on particle production. Azimuthal correlations between 
a forward 7r° and a mid-rapidity particle were first measured with the FPD - and recently 
confirmed, and extended, with the FMS . Comparison of Acp between p+p and d+Au shows a 
significant broadening in the back-to-back peak in d+Au. Such effect appears to be stronger as 
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Figure 2: Left: HIJING impact parameter versus charge sum as recorded by the STAR BBC for simulated 
minimum bias events. Comparison of charge distribution with data is in the inset. Following: uncorrected 
coincidence probability versus azimuthal angle difference between two forward neutral pions in peripheral (center) 
and central d+Au collisions (right). Data are shown with statistical errors and fit with a constant plus two gaussian 
functions (in red). CGC expectations'— ** have been superimposed (in blue) to data for central d+Au collisions. 



the transverse momentum of the particles decreases, as expected by saturation models. One can 
approach the saturation region by lowering the x value of the probed gluon by requiring both 
leading and associated particles to be detected in the forward region, now possible due to the 
FMS wide acceptance. Results from this analysis, where the lowest x is probed, are shown. 

3 Results and Systematics 

Centrality averaged azimuthal correlations between two forward ir° candidates have been mea- 
sured for different px selections. Results are shown on the right-hand side of figure HJ Neutral 
pions are reconstructed from pairs of photon clusters, found within the FMS fiducial volume 
(2.5 < T] < 4.0), that present an invariant mass in the interval 0.05 < M 77 < 0.25 GeV/c 2 . 
The pair with the largest pt is selected as the leading (trigger) 7r° and its azimuthal coordinate 
is compared inclusively with those of all the other (associated) 7r° candidates. The (efficiency 
uncorrected) probability to find an associated 7r° per triggered event presents two contributions. 
The peak centered in Acp = (near-side peak) represents the contribution from pairs of neutral 
pions from the same jet. It is not expected to be affected by gluon saturation, hence providing us 
with a useful tool to check the effective amount of the broadening in the away-side contribution. 
This second peak, centered at Aip = n, represents the back-to-back contribution to the coinci- 
dence probability and it is expected to disappear in going from p+p to d+Au if saturation sets 
in. Data are fit with a constant plus two gaussians centered at Aip = and Aip = ir respectively. 
Figure Q] shows how the width of the near-side peak is not changing from p+p to d+Au, while 
the away-side peak presents significative broadening, with an effect larger than that found in 
forward + mid-rapidity particle correlations^. The pt dependence of the broadening effect is 
studied by applying two different selections to the data. A lower px cut for both trigger (2.0 
GeV/c< 4 r9) ) 



and associated pions (1.0 GeV/c< p^ s < p^ 9 ) shows a stronger broadening 



of the signal width (adAu — &pp = 0.52 ± 0.05) than a more restrictive cut (2.5 GeV/c< p T 



(trg) 



CdAu — &pp = 0.11 ± 0.06), as expected from saturation models. 



1.5 GeV/c< 4 assc) <4 r9 \ 

Broadening effects are expected to be more significant when the more central part of the 
nucleus is probed. In order to disentangle peripheral from central collisions, the sum of charges 
(J2Qbbc) has been recorded using the east side of the STAR Beam-Beam Counter (BBC) 
that faces the Au beam. This provides a measure of the multiplicity of the event, which is 
correlated with the impact parameter in the collision, as shown on the left-hand panel of figure 



[2] for minimum bias events, as simulated by HIJING 1.383^. The multiplicity selections for 
peripheral (0 < J2 Qbbc < 500) and central d+Au collisions (2000 < J2 Qbbc < 4000) result 
in {b} = 6.8 fm and (b) = 2.7 fm respectively. Figure [2] also shows azimuthal correlations for the 
lower px selection, using the charge sum to discriminate between peripheral (central panel) and 
central (right-hand panel) d+Au collision. Peaks on the near-side appear nearly unchanged from 
p+p to d+Au, and particularly from peripheral to central d+Au collisions. On the contrary, the 
away-side peak shows strong differences from peripheral to central d+Au collisions. Peripheral 
d+Au collisions show back-to-back peak like in p+p, even though it appears to be smaller relative 
to the near-side peak than in p+p. Central d+Au collisions show instead large broadening of 
the away-side peak, effectively causing its disappearance. 

Theoretical expectations for azimuthal correlations in interactions between a dilute system 
(deuteron) and a dense target (Gold) have been explored. Inclusive particle production has 
been calculated within the CGC framework, using a fixed saturation scale Qs ■> an d then 
used to compute the di-hadron correlations LLW. Calculations consider valence quarks in the 
deuteron scattering off low-x gluons in the nucleus with impact parameter b = 0. A preliminary 
comparison with data in central d+Au collisions is shown on the right-hand plot of figure [2j 
Calculations (in blue) have been superimposed to the data after adding a constant offset to 
emulate the background from underlying event. CGC calculations show qualitative consistency 
with data in their expectations of a strong suppression of the away-side peak in central d+Au 
collisions. Comparison between same data and CGC calculations using a different approach can 
be found in this reference!^. 

Further systematic studies are being performed for this analysis. Systematic Pythia in- 
vestigation found p+p data to be consistent with gluon distribution function that include a 
rapid rise of the gluon density at low x. Azimuthal correlations were also studied in embedded 
Pythia+GSTAR events into minimum bias d+Au data in order to rule out the possibility that 
additional multiplicity in d+Au events compared to p+p could cause loss of correlation. 

4 Conclusions 

Thanks to a rich d+Au RHIC run in 2008, the Forward Meson Spectrometer is pursuing its 
primary objective of mapping the boundaries for saturation signatures for back-to-back jet 
correlations as a function of r/ and px- A new, interesting piece is provided by correlations 
between two forward neutral pions, which show a strong suppression in the away-side peak in 
central d+Au collision compared to p+p, qualitatively consistent with CGC expectations. 
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